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' Summary 

a • 

The effects which star cluster concentration and binarity have on observable parameters, that 
characterise the dynamical state of a population of stars after their birth aggregate dissolves, 
are investigated. To this end, the correlations between ejection velocity, binary proportion, 
mean system mass, binary orbital period and mass ratio are quantified for simulated aggregates. 
These consist of a few hundred low-mass binary and single stars, and have half-mass radii in the 
range 2.5 to 0.08 pc. The primordial binary-star population has a period distribution similar 
to that observed in Taurus- Auriga for pre-main sequence binaries. The findings presented here 
are useful for interpreting correlations between relative locations and proper motions, binary 
properties and masses of young stellar systems within and surrounding star forming regions, and 
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, of stellar systems escaping from Galactic clusters. 



For the low-concentration binary-rich aggregates, the proportion of binaries decreases 
monotonically as a function of increasing ejection velocity after aggregate dissolution, as 
expected. However, this is not the case for initially highly concentrated binary-rich aggregates. 
The reason for this difference is the interplay between the disruption of binary systems and the 
r/3 • initial depth of the potential well from which the stellar systems escape. 

After aggregate dissolution, a slowly expanding remnant population remains. It can have a 
high binary proportion (80 per cent) with a high mean system mass, or a low binary proportion 
(less than about 20 per cent) with a low mean system mass, if it was born in a low- or a 
high-concentration aggregate, respectively. It follows that adjacent regions on the sky near some 
star-forming clouds can have young populations with different binary proportions and different 
mass functions, even if the binary proportion at birth and the initial mass function (IMF) were 
the same. 

Binary systems that are ejected from the aggregate tend to be massive, and their mass 
ratio tends to be biased towards higher values. The mean system mass is approximately 
independent of ejection velocity between 2 and 30 km/s. Dynamical ejection from binary-rich 
aggregates adds, within 10 Myr, relatively massive systems to regions as far as 300 pc from 
active star- forming centres. Long-period systems cannot survive accelerations to high velocities. 
The present experiments show that a long-period (> 10 4 d) binary system with a large velocity 
(> 30 km/s) cannot be ejected from an aggregate. If such young systems exist, then they will 
have been born in high- velocity clouds. 



Subject headings: stars: formation - stars: pre-main-sequence - stars: low-mass, brown dwarfs - 
stars: kinematics - stars: statistics - binaries: general 
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1. Introduction 

Stellar systems (i.e. single or multiple stars) form in groups. The dynamical processes within these alter 
the properties of the young systems when they leave the site where they formed. The dynamical properties 
of a stellar system are its mass (i.e. luminosity if age is known), the multiplicity and the orbital parameters 
if it is a multiple system. The distribution of velocities of young stars emanating from star-forming centres 
(i.e. the kinematical signature of star formation) will also be affected by the dynamical interactions within 
the young groups. Both, the distribution of dynamical properties and the kinematical signature of star 
formation bear an imprint of the dynamical configuration at the time when the stellar group was born. 

Star formation in Taurus-Auriga gave birth to aggregates with sizes of roughly 0.5-1 pc consisting of 
about 20-50 stars. It is now well established that most stars form in binary systems in Taurus-Auriga (e.g. 
Kohler & Leinert 1998). The same appears to hold true in other star-forming regions (Ghez et al. 1997). 
Embedded clusters may also have a binary proportion that is higher than in the Galactic field (Padgett, 
Strom & Ghez 1997). In the Trapezium cluster, which is a very dense embedded cluster and probably less 
than 1 Myr old, Prosser et al. (1994) find a binary proportion that is at least as large as in the Galactic 
field. In the cluster core, Petr et al. (1998) observe, for low-mass stars, a binary proportion similar to 
the Galactic field, and smaller by about a factor of three than the binary proportion in Taurus-Auriga. 
These findings are particularly interesting, because binary destruction is expected to be efficient in such an 
environment. A review of pre-main sequence binary stars, and their relation to Galactic field systems, is 
provided by Mathieu (1994, see also Kroupa 1995a, Simon et al. 1995). Young Galactic clusters also contain 
binary systems. The particularly well studied Pleiades and Praesepe clusters have binary proportions of 
40-50 per cent, for systems of spectral type earlier than K0 (Raboud & Mermilliod 1998a, 1998b). 

There exists thus evidence that the formation of binary systems may be by far the dominant 
star-formation mode in both loose groups and highly concentrated embedded clusters, some of which 
may evolve to bound Galactic clusters. The term aggregates is used henceforth to mean loose groups or 
embedded clusters of more than 10 stars. 

If stars form predominantly in aggregates of binary systems, then the kinetic energy distribution 
after aggregate dissolution should be enhanced at high energies, when compared to dissolved aggregates 
of single stars, because binary star binding energy can transform into kinetic energy (Heggie 1975, Hills 
1975, Hut 1983). Large accelerations are destructive to binary systems, so that the proportion of binaries 
should be a decreasing function of increasing final kinetic energy. Additionally, different initial aggregate 
concentrations lead to different final binary proportions and kinematical signatures, as will be shown here. 
This is also true under the extreme assumption that all stars always form in binary systems with the same 
initial dynamical properties. This dynamical mechanism of producing variations in binary proportion and 
associated dynamical properties stands in contrast to a possible variation of these parameters determined by 
the star-formation process. Durisen & Sterzik (1994) make the interesting point that the binary proportion 
may be smaller in molecular clouds with a higher temperature than in lower-temperature clouds. 

It is important to study the signatures that arise from purely dynamical interactions in stellar groups, 
for a comparison with outcomes from usually less well-understood alternative scenarios. In a study of the 
large-scale distribution of young stars around active star-forming regions, Sterzik & Durisen (1995) find 
that the dynamical decay of small stellar groups can lead to sufficiently large velocities to populate large 
areas on the sky with young stars, so that these need not have formed near their observed location. They 
find that special initial dynamical configurations of the stars (e.g. cold thin strings) lead to enhanced 
production of ejected stars. Initial decay of cold sub-groups within larger complexes also has this effect 
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(Aarseth & Hills 1972), and scattering of proto-stars on cloud clumps during an even earlier dynamical 
phase may likewise eject very young low- mass stars (Gorti & Bhatt 1996). However, the number of ejected 
stars cannot account for the observed number of widely distributed young stars (Feigelson 1996). The 
evolution of circum-stellar discs around stars ejected from small stellar groups is studied by Armitage & 
Clarke (1997), and McDonald & Clarke (1995) show that the presence of circum-stellar material in small 
proto-stellar groups increases the number of binaries formed and randomises the mass-ratio distribution. 
That the number of dynamically ejected stars is increased significantly in binary-rich stellar aggregates, 
when compared to clusters consisting initially only of single stars, is shown by Kroupa (1995c). These 
simulations show that a mass-ratio distribution produced by randomly associating masses from the IMF, 
decays to the observed distribution for G-dwarf binaries, if most stars form in aggregates similar to observed 
embedded clusters. Also, initially more concentrated aggregates produce more stars with a high ejection 
velocity, the maximum of which increases with decreasing cluster radius. De la Fuente Marcos (1997) 
investigates the dependence on cluster richness, and finds that the mean ejection velocity increases for 
more initially populous clusters. Ejection velocities larger than a few hundred km/s can be achieved in 
young star clusters containing massive primordial binaries (Leonard & Duncan 1990). This may explain 
the location of OB stars far from active star-forming sites. Leonard (1991) finds, on the basis of many 
scattering experiments, that the maximum ejection velocity is of the order of the escape velocity from the 
stellar surface of the most massive star. If its mass is 60 M e , then a similar star can attain an ejection 
velocity of up to 700 km/s. A low-mass star may find itself fleeing with a velocity of up to 1400 km/s, after 
a surface-grazing encounter with such a star. A critical discussion of the possible origin of runaway OB stars 
is provided by Leonard (1995). He stresses that collisions of two stars during binary-binary interactions can 
produce runaway OB stars with very similar properties as in the alternative scenario, in which such stars 
result from a supernova explosion in close binary systems. An interesting and insightful discussion of the 
implications of the binary properties of runaway OB stars for the dynamical configuration of massive stars 
at birth is to be found in Clarke & Pringle (1992). 

In this paper, the correlations between stellar velocity, system mass and binary proportion that arise 
from aggregates with different initial concentration and consisting initially either of 400 single stars or of 
200 binary systems, is studied. The resulting correlations are useful for interpreting the properties and 
distribution of young stars near and in star forming regions (see for example Brandner et al. 1996, Feigelson 
1996, Frink et al. 1997). 

In Section || the assumptions, simulations and definitions are described. The results are presented in 
Section 0, and Section contains the conclusions. 



Method 



The initial conditions and numerical method are described in Section 2A, and the data analysis is outlined 
in Section 2.2. 



2.1. Assumptions 

Abin = 200 binary systems are distributed in virial equilibrium according to the Plummer density law, with 
initial half mass radii i?o.5 = 2.53,0.77,0.25,0.077 pc. These approximately span the region of parameter 
space similar to distributed (e.g. Taurus-Auriga) and very tightly clustered (e.g. Trapezium cluster) 
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star formation. The clusters have zero initial ccntre-of-mass velocities in the local standard of rest. The 
Ro.5 — 0.8 pc aggregate is especially interesting, because inverse dynamical population synthesis (Kroupa 
1995a, b) suggests that it may be representative of the dynamical structures in which most stars form 
(compare with Lada & Lada 1991). 

While the mechanism of binary system formation cannot be specified in detail, the assumption that 
the majority of all stars form in binaries is supported by observational evidence (see review by Mathieu 
1994), and by recent advances in the theory of star formation (for reviews see Boss 1995, Clarke 1996). 
However, theory cannot, at present, constrain the early dynamical properties of stellar systems. The 
interesting suggestion has been made (Durisen & Sterzik 1994) that cloud temperature may influence the 
binary proportion, such that it may be lower in dense embedded clusters. For comparison with the binary 
rich aggregates, iV s i ng = 400 single stars are distributed in aggregates, initially with i? .5 = 0.25,0.077 pc. 

The initial velocity dispersion, a, and escape velocity from the centre of the aggregates, w osc = y/2 \<j>\ (</> 
is the Plummer potential at the origin), are: a — 0.3 km.s, v csc — 0.77 km/s (R0.5 = 2.5 pc), a = 0.5 km/s, 
w osc = 1.4 km/s (R0.5 — 0.8 pc), a = 0.9 km/s, v csc = 2.4 km/s (R0.5 = 0.25 pc), and a = 1.7 km/s, 
w osc = 4.4 km/s (Rq.5 = 0.08 pc). Other physical parameters are listed in table 1 of Kroupa (1995a). 
Aarseth's NBODY5 programme (Aarseth 1994) is employed for the N-body simulation of the dynamical 
evolution of each aggregate in a standard Galactic tidal field. 

In order to simplify the computational burden, the stars are treated as point particles and stellar 
evolution is neglected. The assumption of virial equilibrium is the simplest case, and implies that the 
results presented here are strictly only applicable to escaping stars from Galactic clusters. The present 
results can, however, also be used as guidelines of the type of correlations one might find after embedded 
clusters dissolve. An explicit formulation of this problem requires treatment of gas expulsion, and thus 
the introduction of additional ill-defined parameters. As gas expulsion is not treated here, the results are 
representative of star formation with high efficiency, i.e. aggregates with low residual gas content. In the 
alternative case of a low star formation efficiency, the major effect gas expulsion has, is a shortening of the 
time-scale during which the dynamical evolution occurs. This can be compensated for by a reduction of 
R0.5, in order to obtain the same effective dynamics (section 6.4 in Kroupa 1995a). The initial «esc is then 
larger. 

Stellar masses, m, with 0.1 Mq < m < 1.1 Mq, are obtained from the IMF: £(m) oc mT a \ a.\ = 1.3 for 
0.08 M Q < to < 0.5 Mq, a 2 = 2.2 for 0.5 M© < to < 1.0 M© (Kroupa, Tout & Gilmore 1993), and a 3 = 2.7 
for 1.0 Mq < m (Scalo 1986), where £(m) dm is the number of stars with masses in the range to to m + dm. 
The mean stellar mass is 0.32 Mq, and the mass of each aggregate amounts to M tot = 128 Mq. Adopting 
for the mass of B stars 6 — 18 Mq, each should have associated with it 280 stars with mass in the range 
0.08 — 1 Mq. The maximum ejection velocity that can be achieved is thus limited to about 600 km/s for 
0.1 Mq stars, and about 300 km/s for G-dwarfs (Leonard 1991). 

The main-sequence mass-ratio distribution for G-dwarf binaries (Duquennoy & Mayor 1991) is not 
consistent with random pairing from the IMF, but may be derived from this assumption if most stars form 
in embedded clusters (Kroupa 1995a, b). In accordance with this result, and the evidence presented by 
Leinert et al. (1993), stellar masses are combined at random to generate the initial binary-star population. 
Special care must be taken when interpreting an observed mass-ratio distribution, as it can be affected 
significantly by even simple observational bias (Trimble 1990, Tout 1991). 

Binary systems must arrive on the birth-line with eccentricities approximately dynamically relaxed, 
because subsequent thcrmalisation in the stellar aggregate is not efficient enough to produce such a 
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distribution (Kroupa 1995b). This is because the cross-section for a significant change in eccentricity 
decreases very steeply with increasing distance of closest approach of a perturber (Heggie & Rasio 1996). 
Consequently, the initial eccentricity distribution is taken to be dynamically relaxed. The results are not 
sensitive to this assumption, however. 

An initial period distribution that is consistent with the observational data for young binaries is used. 
The orbital periods, P (in days), form a flat distribution, /p(log 10 P) = [log 10 (P max ) — log 10 (P m i n )] _1 
(equation 3 in Kroupa 1995a), with logio-P m in = 3, logi P m ax = 7.5 and P m { n < P < P ma x- 

For each binary and single-star aggregate, A run = 5 and 3 simulations, respectively, are carried through. 



2.2. The observables 

All results quoted here are averages of A mn simulations that are evaluated after 1 Gyr, i.e. after the 
aggregates have dissolved. Aggregate dissolution occurs after 700 ±130 Myr in all cases, when the number 
of stars in a volume with a radius of 2 pc, that is centred on the density maximum of the cluster, has 
decayed to 3 of less. 

The velocity-dependent binary proportion is 

t - Nhin ' v m 

Jv (N- +Nw V {) 

where A s i ngjt , and Abi n ,u are the number of single-star and binary systems, respectively, in a velocity interval 
v to v + Av relative to the local standard of rest. The binary proportion in some sub-domain, which may, for 
example, be the period range or spatial region accessible to the observer, is / = Abi n /(A s i ng + Abi n ), wnere 
A s i n g and Ab; n are the number of single and binary systems, respectively, in the sub-domain. Similarly, / to t 
is the binary proportion of the entire population. 

The mean system mass in the velocity interval is 

<m> v=m — Tn — v (2) 

l/ v sing,i> ~T JVbm,v ) 

where M v is the total stellar mass in the velocity interval. Initially, i.e. at t = 0, /„ = 1 and 

< m >„= 0.64 Mq independent of v for the binary-star aggregates, and f v = with < m > v = 0.32 M Q for 

the single-star aggregates. 

The relative proportion of systems in a velocity interval is 

^ _ ( Aging. t, + Abi n ,t>) ^ 
(Aging, tot + Abin,tot) 

Note that this is f v in Kroupa (1995c). 

The circular orbital velocity, u or b [km/s], of a binary star with system mass, m sys [Mq], and orbital 
period, P [days], is 
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log 10 P = 6.986 + log 10 m sys - 3 log 10 u or b. (4) 

The primordial binary population used here has a maximum u or t, = 27.7 km/s (for P — 10 3 d, 
TO sys = 2.2 Mq) and a minimum v OT b = 0.39 km/s (P = 10 7 5 d, m sys — 0.2 M Q ). 

Finally, a stellar system that is ejected has a final (i.e. evaluated after 1 Gyr) velocity v > 2 km/s. 

3. Results 

Long distance encounters between systems (two-body relaxation) and in addition scattering of systems 
on the non-uniform background potential (collective effects) dominates the dynamical evolution of the 
aggregates. Relatively energetic, stochastically occurring encounters between stellar systems can lead to the 
ionization of binary stars and to the acceleration of a system to escape velocity from the aggregate. 

If the mean kinetic energy of the population is -Ekin and the binding energy of a binary is 
— -Efcin = —Gmim2/(2a), where m t and a are the masses of the components and the semi-major axis, 
respectively, then a binary is termed hard if Ebin/E^n > 1- Hard binaries are likely to gain binding energy 
i.e. to harden (Heggie 1975, Hills 1975), in which case the perturber can be accelerated to escape velocity. 
The resulting hardened binary suffers a recoil which may be sufficient to also expel it from the aggregate. 

These processes change the distributions with velocity of the number of systems, h v , of the binary 
star proportion, f v , of the mean system mass, < m > v . Also, the correlations between binary-star binding 
energy, E^ n , and kinetic energy, -Ekim an d between the velocity, v, and orbital period, system mass and 
mass ratio, evolve. The distributions that emerge after aggregate dissolution thus contain information 
about the initial dynamical configuration, but care must be taken in interpreting distribution data. This is 
the subject of the present section. 

3.1. Distribution of velocities 

In Fig. [j], h v , f v and < m > v are plotted as a function of velocity for the binary star aggregates with 
Ro.5 = 2.5,0.8,0.25,0.08 pc. Fig. || contains the same information for the two single star aggregates. In 
Table |l|, column 1 contains the centre of each logarithmic velocity bin. Columns 2-7 list, for each aggregate, 
the fraction, h v , of systems per logarithmic velocity bin. 

After aggregate dissolution, most systems have a velocity near 0.35 km/s (Figs. @ and §). A slight 
shift in the maximum of h v towards smaller v with decreasing i?o.5 comes about, because systems have 
to overcome the initial aggregate potential before escaping. The cooling is much more pronounced in the 
absence of binary star heating (Fig. ||), and during the first few cluster crossing times. Later, the aggregate 
expands to fill its tidal radius and looses memory of its initial concentration. Initially very concentrated 
aggregates, and those with large R0.5, have an indistinguishable life-time (Kroupa 1995c). The binary star 
population, however, retains this memory (Kroupa 1995a). 

For aggregates with initially smaller -R0.5, an increase in the proportion of systems with v > 2 km/s 
results. Comparison of Figs. |] and || shows that, for the same i?o.5, aggregates initially with a high 
proportion of binary systems have significantly larger h v at v > 2 km/s, than aggregates that consist 
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initially only of single stars. A high proportion of primordial binary systems thus increases the percentage 
of ejected systems. 

3.2. Binary proportion 

High velocity systems are expected to be primarily single stars, because only relatively hard binaries can 
survive the large accelerations during the encounter, as is also stressed by Sterzik & Durisen (1995). This 
is borne out for the i?o.5 = 2.5,0.8,0.25 pc aggregates (Fig. [l]), and also in the simulations reported by 
Leonard & Duncan (1990, their figs. 4 and 5) and also fig. 9 in Kroupa (1995c). 

Dynamical evolution is quiescent in star forming regions where the stellar systems freeze out of 
the gas in low-density aggregates. The i?o.5 = 2.5 pc aggregate approximates this situation. In this 
aggregate, a binary with m sys = 0.2 M Q and P = 10 7 5 d has v Q tb — 0.39 km/s, which is comparable to 
the velocity dispersion. The entire binary population is therefore hard, and indeed most binaries survive 
cluster evolution. From such aggregates results a high (f v > 0.8) proportion of binaries for systems with 
approximately v < 1 km/s. Only 2.6 per cent of the systems end up with v > 5 km/s (Table |l|), and these 
have a binary proportion of approximately f v < 0.1 (Fig. |l|). 

For the i?o.5 = 0.8 pc aggregate, a relatively high proportion of binaries (/„ > 0.8) among systems that 
have v < 0.3 km/s is obtained. Significant differences between the binary proportions in the two models 
can be found in the velocity intervals a) —0.5 < logi w < and b) 0.1 < logi w < 0.6. In interval (a), 
f v ~ 0.9,0.6 and in interval (b), f v k 0.05,0.15, for i? .5 = 2.5,0.8 pc, respectively. Of all systems that 
finally emerge from such an aggregate, 4 per cent have a velocity v > 5 km/s. These have a slightly larger 
binary proportion than for the i?o.5 = 2.5 pc case discussed above. Overall, for i?o.5 > 0.8 pc, f v decreases 
monotonically with increasing v, and escaping stars have a low (/„ < 0.2) binary proportion. 

Star formation in denser aggregates (i?o.5 < 0.8 pc) leads to significantly different behaviour of /„ with 
v (Fig. [I]). For i? .5 = 0.25 pc, f v « 0.45 for v < 1 km/s, with a discontinuous decrease to f v < 0.2 for 
v > 1 km/s. 

Of special interest is the i?o.5 = 0.08 pc model. It represents most closely the Trapezium cluster, 
because it has a comparable central number density, half-mass radius and velocity dispersion. In the present 
model, the initial crossing and relaxation times are t cx = 1 x 10 5 yr and t le i — 3 x 10 5 yr, respectively 
(Table 1 in Kroupa 1995a), whereas in the Trapezium cluster, i cr k 4 — 12 x 10 5 yr and t TC \ « 0.7 — 3.7 Myr 
(Bonnell & Kroupa 1998). The Trapezium cluster, however, is different in that it contains 500-1000 stars 
with a mean mass of about 0.6 M@, and in that it is the core of the much more massive and extended Orion 
Nebula Cluster (Hillenbrand & Hartmann 1998). Also, it is not clear if the entire cluster is gravitationally 
bound. 

After dissolution of this R0.5 — 0.08 pc aggregate most systems have v < 1 km/s (Fig. |l|). The 
proportion of binary stars shows a rather complex dependence on v. The binary proportion ranges from 
/„ rs 0.1 for systems with v ~ 0.1 — 0.2 km/s to /„ ~ 0.6 for v 1 km/s; /„ thus increases with v 
for v < 1 km/s. The binary proportion shows a significant maximum (f v w 0.6) near v m 1 km/s, and 
remains approximately constant at f v ss 0.25 for 3 km/s< v < 20 km/s. In this rather extreme model of 
star formation, 4.6 per cent of all systems have v > 5 km/s after aggregate dissolution. The low value 
of f v for small v, and its rise with v, is due to efficient disruption of binaries at an early dynamical age, 
when the ejected stars are decelerated most effectively by the young deep potential well. A fraction of 



- 8 - 



the predominantly single stars with low velocity, is also a decelerated part of the high-velocity tail in the 
aggregates with R > 0.8 pc. This can be inferred from Figs. || and ^ (Section 3.4). Hardened binaries 
(logio-P < 4) are found with small velocities. Usually they are the result of energetic three-body or 
binary-binary interactions causing ejection. 

Thus, in a realistic embedded cluster with the same stellar mass, the initial escape velocity is larger 
because the potential is dominated by the gas. Within a few Myr most of the gas is expelled, leaving an 
expanding cluster population, and a binary deficient remnant population, in which each system has a small 
centre-of-mass velocity. This decelerated and binary deficient population remains bound to the molecular 
cloud and, after a few Myr, contributes to a distributed population of young stars with significantly different 
dynamical properties to the distributed population in Taurus-Auriga. Dispersal of this binary deficient 
population takes long, and an observer finds a loosely distributed group of stars of similar age, and with a 
reduced binary proportion that depends on the initial cluster concentration. 

Concerning the aggregates with initially no primordial binaries (Fig. |^) , more binaries form by capture 
in the initially more concentrated aggregate (i?o.5 = 0.08 pc), owing to the more frequent three-body 
encounters in the young concentrated aggregate. The resulting total binary proportion, however, remains 
insignificant (Kroupa 1995a). The data plotted in Fig. || indicate that /„ increases with v (for v > 1 km/s), 
which is contrary behaviour to the aggregates that contain a large population of primordial binaries. 



3.3. Mean system mass 

In a binary-binary or binary-single star encounter, binding energy can be transformed into kinetic energy 
of the escaping stellar system. A given acquired kinetic energy corresponds to a smaller ejection velocity if 
the system mass is larger. Low-mass stars can be ejected with higher velocities than high-mass stars. The 
extensive simulations performed by Leonard & Duncan (1990) and Leonard (1991, 1995) demonstrate that 
this is the case, and are consistent with the observational mass-velocity diagram for OB stars produced by 
Gies & Bolton (1986). This is also discussed in length by Conlon et al. (1990). 

The present study concentrates on the mean-system-mass-velocity relationship obtained from 
self-consistent A-body simulations of clusters of low-mass stars (m < 1.1 M Q ), and is thus relevant for 
the large-scale distribution of young-low mass stars seen in the ROSAT survey (compare with Sterzik & 
Durisen 1995). 

As is evident from Figs. [I] and ^, the behaviour of < m > v with v is complex and depends on the initial 
concentration of the aggregate. The complexities of the underlying binary-binary and triple-star encounters 
are discussed at length by Harrington (1975), Heggie (1975), Hills (1975), Leonard & Duncan (1990) and 
Leonard (1991, 1995), and a review can be found in Valtonen & Mikkola (1991). The most-apparent result 
here is that the expected simple correlation (smaller < m > v for larger v) does not hold, except for the 
i?o.5 = 2.5 pc aggregate. The few stars that are ejected from this aggregate are low-mass stars expelled 
from unstable triple or higher-order systems (Heggie 1975, Harrington 1975, Hills 1977, compare also with 
the second-mass-family interactions of Leonard 1991, and with Kiseleva, Eggleton & Orlov 1994). 

The expected correlation is also observed for v > 30 km/s for all binary-star aggregates. An acceleration 
beyond this velocity, which is the orbital velocity of the hardest primordial binary in the present simulations 
(Section |2.2| ), is destructive to all primordial binaries. Only single stars appear with such large velocities. 

For i?o.5 = 0.8,0.25 and 0.08 pc, the velocity range v w 2 km/s to 30 km/s yields no clear correlation 
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between < m > v and v: < m > u « 0.5 Mq is approximately constant. This is an interesting finding which 
was also noted by Harrington (1975). It shows that systems more massive than 0.5 Mq are ejected from 
aggregates similar to embedded clusters, with velocities that can place them at distances between 20 and 
300 pc within 10 Myr of ejection time from their formation site. These systems result from stochastic and 
quite energetic binary-binary and three-body encounters. 

The majority of systems that have v < 1 km/s, and which do not spread much further than 1-10 pc 
from their formation site within 10 Myr, have a constant < m > v that lies between the average stellar 
mass and mean primordial system mass. Smaller values are seen for binary-poor remnant populations 
(e.g. for i?o.5 = 0.08 pc). It follows that an observer must be careful not to interpret the stellar mass 
function and binary proportion of such populations in terms of a possible dependence of these quantities on 
star- formation environment, without due consideration of the dynamical history. 

3.4. Binding energy, period and mass ratio 

Only binaries that are sufficiently bound will not be ionised when they suffer a close interaction with 
another system, after which they may leave the aggregate with relatively high velocity. Thus, the correlation 
between binary star binding energy, -Ebin, and its centre of mass kinetic energy, -Ekin, indicates the history 
of a system. 

Above it was seen that the final binary proportion, /„, for systems with 2 km/s< v < 30 km/s, 
is larger for initially more concentrated binary star aggregates. It can achieve values of 20-40 per cent 
for i?o.5 < 0.25 pc, although the overall final binary proportion is small (/tot ~ 0.27, fig. 3 in Kroupa 
1995a). This result is relevant for the distribution of young stars around star-forming regions. Initially 
highly concentrated embedded clusters may add young binaries to regions as far as 300 pc over a period of 
10 Myr. Such binaries have a well-defined correlation between E^i n an d -^kin- This correlation transforms 
to correlations between P and v, between the mass ratio (q = 777,2/mi < 1) and v, and between the system 
mass (mi + 777.2) and v, where mi and 7712 are the primary- and secondary-star masses, respectively. 

3.4-1- Binding energy and period 

In Figs. H and [| are plotted i?bin against Pkin, as well as the orbital period, P, against ejection velocity, 
v, for each binary system in the -/V run simulations. The distribution of data points for the i?o.5 = 2.5 pc 
aggregate (top two panels in Fig. ^) reflects approximately the initial distribution in binding energies 
(logio-Ebin > 2) and orbital periods (logio-P > 3). Only very few binaries have gained binding energy and/or 
have been accelerated to higher velocities. As the initial R0.5 is reduced, the number of binary systems 
with hardened orbits and larger v increases. The depletion of orbits at large P is clearly evident in the 
R0.5 — 0.08 pc cluster. 

Figs. [| and ^ nicely show that a binary with an orbital period corresponding to t; or b only remains bound 
when suffering a collision, if the ejection velocity v < v OI b. A single notable exception, that occurred in the 
five simulations of the R0.5 = 0.08 pc aggregate, is the binary system with P = 10 5 7 d and v = 16 km/s. 
It is difficult to trace the detailed dynamical history of any individual stellar system owing to the discrete 
output times when stellar masses, positions and velocities are written to computer disk. But this binary 
probably formed in a complex high-order interaction, that resulted in two stars being ejected on essentially 
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identical trajectories. 

The results discussed so far are valid for a primordial binary star population that has periods P > 10 3 d. 
In reality, binaries with shorter periods do exist with a binary proportion / « 0.15 (fig. 1 in Mathieu 1994). 
The inclusion of primordial binaries with P < 10 3 d does not change the results presented here, apart from 
slightly increasing f v for the ejected systems. In Figs. || and ||, both the (.Ebin, E^ n ) and (logioP, v) plots 
would contain orbits with i?bin > 10 2 M km 2 /s 2 and P < 10 3 d, and an upper envelope for v > 1 km/s 
given by the dashed diagonal lines. 

The correlation between binding and kinetic energy is established particularly well for the binary 
systems that are formed by capture in the single star aggregates (Fig. ||). Their periods range from 
P > 100 d to 10 11 d. The binary systems with P m 10 10 d form by chance three-body low-velocity 
encounters and form a distinct group in the figure. 

3.4-2. Mass ratio and system mass 

More massive systems have a higher binding energy, and are thus more resistant to higher accelerations. 
Binaries that are ejected from an aggregate should thus have a larger system mass and a mass-ratio nearer 
to unity for higher v. 

In Figs. ^ and mi + m% and q are plotted in dependence of v for the binary-star aggregates. The 
Ro.5 = 2.5 pc aggregate shows approximately the initial distributions, and is useful as a reference. As the 
aggregate concentration is reduced, the number of ejected systems increases. These tend to have larger 
system masses, as expected. At the same time, binary stars with a mass-ratio q < 0.2 are preferentially 
removed as i?o.s is reduced. This is expected because, for a given semi-major axis, they have the lowest 
binding energy, -Ebin oc m 2 x q. The correlation between q and v is in the expected sense, which is evident 
in the figures by the appearance of orbits in the region q > 0.6, v > 2 km/s. However, the correlation 
is weak, because -Ebin also depends on mi. As is evident from the figures, the number of binaries with 
mi +77i2 > 1-5 Mq increases for smaller i?o.5- This is because the more frequent three-body and higher-order 
interactions lead to more frequent exchanges of companions, which most often leads to the production of 
binary systems consisting of the most massive stars involved in the interaction (Harrington 1975, Heggie 
1975, Hills 1977, see also Valtonen & Mikkola 1991, McDonald & Clarke 1993). 

Concerning the single-star aggregates, Fig. || shows that the distribution of mass ratios of the 
dynamically formed binaries in the i?o.5 — 0.25 pc aggregate is roughly flat over the whole accessible range. 
In the more concentrated aggregate, binaries form with larger system mass and a bias towards larger q, 
which is a result of dynamical biasing discussed in greater detail by McDonald & Clarke (1993). 

4. Conclusions 

The correlations between ejection velocity and the proportion of binaries, as well as their orbital parameters, 
have been quantified for a range of initial dynamical configurations. The correlations are useful in the study 
of stellar systems that are apparently ejected from Galactic clusters (see e.g. Frink et al. 1997), some of 
which are known to be rich in binaries (e.g. Raboud & Mermilliod 1998a, 1998b). Observed ejected binaries 
should show correlations as presented in Figs. [| Ui §| an d @- The results for the binary-rich aggregates 
modelled here are also relevant for an understanding of the large-scale distribution of young stars, because 
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most stars appear to form in aggregates with a high binary proportion. Additionally, the correlations 
contain information about the dynamical configuration at birth. 

For binary-rich aggregates containing a few hundred stars the following correlations result: (i) more 
tightly clustered aggregates lead to more stellar systems having larger ejection velocities and a smaller 
overall binary proportion, (ii) the large population of primordial binaries leads to a significantly enhanced 
number of systems with high-ejection velocities compared to single-star aggregates, (iii) systems with 
high ejection velocities have a significantly reduced binary proportion, (iv) binary stars with high ejection 
velocities have short-period orbits, and tend to be more massive with a mass-ratio biased towards unity, 
(v) the average system mass as a function of ejection velocity is defined above about 2 km/s by stochastic 
close encounters, so that systems more massive than 0.5 Mq with high ejection velocities occur, and (vi) 
aggregates with i?o.5 < 0.25 pc lead to a complex dependence of the resulting binary proportion on velocity, 
whereas a stellar population emerging from less concentrated aggregates shows a monotonic decrease of the 
binary proportion with increasing velocity. 

For aggregates of a few hundred single stars one obtains: (i) more tightly clustered aggregates lead 
to an increased number of stellar systems with larger ejection velocities (but significantly less so than in 
the binary-rich aggregates), and an enhanced overall binary proportion that remains significantly below the 
observed binary proportion in the Galactic field, (ii) the binary proportion increases with ejection velocity, 
(iii) is as (iv) above, and (iv) is as (v) above. 

Remnant unbound young populations take long to disperse because they have a small velocity 
dispersion. The binary proportion and mean system mass (and thus the inferred IMF) of such a remnant 
population, sensitively depends on the initial dynamical configuration of the binary-rich birth aggregate. 
After emerging from the birth aggregate, the distribution of velocities of a young stellar population changes 
with time in the gravitational potential of the nearby molecular clouds. A substantial proportion of 
emerging stars is likely to remain bound to the parent molecular cloud until it ceases to exist. 

These findings are important for interpreting the spatial distribution, kinematics and binarity of young 
stars within and surrounding star-forming regions. Molecular clouds, in which stars form preferentially in 
dense embedded binary-rich clusters, should have an enhanced halo population of ejected and relatively 
binary poor (/ w 0.25) young stellar systems. Also, young but binary-depleted groups of stars can be 
misinterpreted to be evidence for an environmental dependency of the binary-formation mechanism. For 
example, in fig. 6 of Brandner et al. (1996), the region US-B has more binaries than the region US- A, which 
also contains many more B stars than US-B. The presence of B stars suggests that the stars in US- A may 
have formed in dense embedded clusters. The stars seen in US-A would then constitute the v < 1 km/s 
remnant population for i? .5 £ 0.25 pc (Fig. 0). Given the results of the present study, it is suggested 
that such a difference in binary proportion between two regions may be due to different initial dynamical 
configurations, and need not imply a dependence of the binary proportion on the star- forming environment. 

Important for the interpretation of the large-scale distribution of young stars surrounding star forming 
sites is the realisation that relatively massive systems are ejected with relatively large velocity (2-30 km/s, 
Fig. g), which is a point also stressed by Sterzik & Durisen (1995). The X-ray surveys are flux limited and 
detect the massive stars (Wichmann et al. 1996), the presence of which around star- forming regions may be 
a natural consequence of the processes studied here. However, if some young binary systems are found to 
have orbital periods that place them above the dashed lines in the right panels of Figs. then this would 
support the suggestion by Feigelson (1996), that some star-formation occurs in small high-velocity clouds. 
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Table 1: Velocity distributions. R% 5 is for the two single-star aggregates. 
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Fig. 1. — Binary proportion, f Vl as a function of velocity, v, for the aggregates with i? .5 = 2.5,0.77,0.25 
and 0.077 pc (solid histogram). Solid circles are the mean system mass, < m > v . Error-bars are Poisson 
uncertainties. Long-dashed curves represent the distribution of velocities, h v , plotted in arbitrary ordinate 
units and tabulated in Table 1. Horizontal dashed lines mark the rough one-sigma range of the binary 
proportion for G-, K-, and M-main sequence systems in the Galactic field (/ t °^ s = 0.47 ± 0.05, Kroupa 
1995a). From the value of f v around the maximum of h v , as well as the shape of the period and mass- 
ratio distributions, it is deduced that 0.25 pc < i?o.5 < 0.8 pc are solutions to inverse dynamical population 
synthesis (Kroupa 1995a). 
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Fig. 2. — As Fig. |l| but for stellar aggregates that initially contain no binary system. 
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Fig. 3. — Distribution of orbits in the binding-energy-kinetic-energy plane (left panels) (energy units are Mq 
km 2 /s 2 ) and the period-velocity plane (right panels) for the two binary-star aggregates with i?o.5 = 2.5 pc 
and 0.77 pc. Orbits along the diagonal dashed line have a binding energy that is equal to the kinetic energy. 
The region between the two vertical dashed lines for R0.5 — 2.5 pc and R0.5 = 0.08 pc (Fig. ||) contains 95 per 
cent of all initial kinetic energies. The vertical dotted line plotted for these two aggregates approximates 
-Skin- Orbits with E^/E^ > 1 are termed hard and are less likely to be ionised. Dashed lines in the right 
panels are equation 4 for m sys = 2.2 Mq (upper line) and m sys = 0.2 M Q (lower line). The vertical dotted 
lines are the initial velocity dispersions. 
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Fig. 4. — As Fig. || but for the two binary-star aggregates with i?o.5 = 0.25 pc and 0.08 pc. 
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Fig. 5. — As Fig. || but for the two single-star aggregates with i?o.5 = 0.25 pc and 0.08 pc. 
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Fig. 6. — Distribution of orbits in the system-mass-velocity plane (left panels) and the mass-ratio-velocity 
plane (right panels) for the two binary-star aggregates with R0.5 = 2.5 pc and 0.77 pc. The system mass 
is mi + to 2 and the mass ratio is q — m 2 /mi < 1, where mi and m 2 are the primary- and secondary star 
masses, respectively. The vertical dotted lines indicate the initial velocity dispersion. 
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Fig. 7. — As Fig. || but for the two binary-star aggregates with i?o.5 = 0.25 pc and 0.08 pc. 
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Fig. 8. — As Fig. || but for the two single-star aggregates with i?o.5 = 0.25 pc and 0.08 pc. 



